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ABSTRACT 

We estimate the distribution of intrinsic shapes of APM galaxy clusters from the 
distribution of their apparent shapes. We measure the projected cluster ellipticities 
using two alternative methods. The first method is based on moments of the discrete 
galaxy distribution while the second is based on moments of the smoothed galaxy 
distribution. We study the performance of both methods using Monte Carlo cluster 
simulations covering the range of APM cluster distances and including a random dis- 
tribution of background galaxies. We find that the first method suffers from severe 
systematic biases, whereas the second is more reliable. After excluding clusters dom- 
inated by substructure and quantifying the systematic biases in our estimated shape 
parameters, we recover a corrected distribution of projected ellipticities. We use the 
non-parametric kernel method to estimate the smooth apparent ellipticity distribu- 
tion, and numerically invert a set of integral equations to recover the corresponding 
distribution of intrinsic ellipticities under the assumption that the clusters are either 
oblate or prolate spheroids. The prolate spheroidal model fits the APM cluster data 
best. 

Key words: galaxies: clustering - Cosmology: Large Scale Structure in the Universe 



1 INTRODUCTION 

Clusters of galaxies are the largest gravitationally collapsed 
objects in the universe and their internal dynamics and mor- 
phologies provide useful cosmological information. In recent 
years many studies of cluster shapes and orientations have 
showed that they are strongly elongated, maybe more so 
than elliptical galaxies, and they tend to point towards their 
neighbours (Carter & Metcalfe 1980; BingeUi 1982; Di Fazio 
& Flin 1988; Plionis, Barrow & Frenk 1991; De Theije, Kat- 
gert & van Kampen 1995). Plionis, Barrow and Frenk (1991) 
(hereafter PBF) have computed ellipticities and major axis 
orientations for the largest up to date sample of about 400 
Abell clusters and found that their apparent shapes are con- 
sistent with those expected from a population of prolate 
spheroids. Support to the prolate spheroidal case was pre- 
sented recently by Cooray (1999) analysing a sample of 25 
Einstein X-ray clusters of Mohr et al. (1995). Struble and 
Ftaclas (1994) analysed a compilation of 344 Abell clus- 
ter ellipticities and found that rich clusters are intrinsically 
more spherical than poorer clusters. In the same framework 
McMilan et. al (1989) studied the ellipticities and orien- 
tations of 49 Abell clusters using Einstein X-ray data to 
trace the hot gas, and also found that the cluster poten- 



tial is quite flat although less so than that found in optical 
studies. Buote & Canizares (1996) analyzed ROSAT PSPC 
images for 4 Abell clusters (including Coma) and, assum- 
ing hydrostatic equilibrium, they found ellipticities of order 
emaas — 0.40 — 0.55 (see also Canizares & Buote 1997). 

Theoretical expectations regarding cluster shape and 
morphology have been investigated via N-body simulations, 
which show that the intrinsic shapes of simulated clusters 
are rather triaxial with an almost uniform distribution of 
shapes between prolate and oblate spheroids (cf. Frenk et 
al. 1988; Efstathiou et al. 1988). Detailed analysis of clus- 
ter morphological parameters and substructure, utilising the 
concept of power ratios (cf. Boute & Tsai 1994) , can be used 
to constrain different cosmological models (cf. Thomas et al 
1998; Valdarnini, Ghizzardi & Bonometto 1999). 

It is obvious that information about the intrinsic shape 
of a cluster is lost when projected on the plane of the sky. 
Many different studies have attempted to recover the dis- 
tribution of intrinsic cluster shapes from the corresponding 
apparent distribution using inversion techniques based on 
the assumption that their orientations are random. 

In this paper we use the APM catalogue (Dalton et 
al. 1997) to measure the projected shape distribution, cor- 
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rected for various systematic effects, and hence attempt to 
estimate the intrinsic shape of clusters. The APM clusters 
are typically as rich as Abell R = clusters, but due to the 
careful identification procedure do not suffer from significant 
projection effects. 

The plan of the paper is the following: In Section 2 we 
describe the APM galaxy and cluster survey, in Section 3 we 
present our projected cluster shape determination method 
and by using Monte Carlo simulations we establish its statis- 
tical robustness. Wc discuss how the foreground/background 
contamination (projection effects) affects the projected clus- 
ter shapes and present a statistical ellipticity correction pro- 
cedure. The correction assumes that clusters are in dynami- 
cal equilibrium and therefore wc exclude clusters with strong 
substructure. In Section 4 we invert the systematic bias- 
corrected projected ellipticity distribution to recover the in- 
trinsic one. Our conclusions are presented in Section 5. 



2 THE APM DATA 

The APM survey covers an area of 4300 square degrees in 
the southern sky (b < —40°) and contains about 2.5 mil- 
lion galaxies brighter than a magnitude limit of bj = 20.5. 
Details of the APM data can be found in Maddox et al. 
(1990a and 1990b), and Maddox, Efstathiou & Sutherland 
(1996). Here we present only a brief summary of the cata- 
logue. The survey was compiled from 185 survey plates from 
the UK Schmidt telescope scanned by the Automatic Plate 
Measuring (APM) machine in Cambridge. The scanned re- 
gion of each plate covers 5.8° x 5.8° of the sky, and since 
neighbouring plate centers are separated by 5° this leads 
to 0.8° overlaps along plate boundaries. The data for each 
plate is stored separately to preserve the multiple measure- 
ments in the overlap regions. Extensive internal checks and 
external calibration have shown that the plate-to-platc zero- 
point error has an rms of 0.06 magnitudes, and that large- 
scale photometric gradients are even smaller. The IRAS and 
COBE all-sky maps show that the galactic obscuration in 
this region of the sky is typically 0.06 magnitudes introduc- 
ing comparable uncertainty in the photometry. The image 
profiles and shapes were used to classify them into gala^x- 
ies, stars and blended stars. Visual checks and deeper CCD 
images show that the classification leads to galaxy samples 
which are 90-95% complete with contamination of 5-10% 
from non-galaxies. 

Dalton et al (1997) applied an object cluster finding al- 
gorithm to the APM galaxy data, and so produced a list 
of galaxy clusters, most of which have subsequently been 
spectroscopically confirmed as clusters. The cluster finding 
algorithm consists of two main steps: The first step uses a 
percolation algorithm to link all pairs of galaxies with sep- 
arations < 0.7 the mean inter-galaxy separation. All mu- 
tually linked pairs arc joined together to form groups, and 
the groups with more than 20 galaxies are identified as can- 
didate clusters. In the second step, an iterative routine is 
applied to each candidate to estimate the richness and char- 
acteristic apparent magnitude of galaxies within a search 
radius of 0.75 Mpc. This produced a list of 957 clusters 
with Zest'^ 0.1 and APM richness of more than 40 galaxies, 
corresponding roughly to Abell richness class 0. The angular 
diameter of the search radius is set to be consistent with the 



distance estimated from the apparent magnitude of galaxies 
in the search radius. 

For our present analysis we cross-correlated the cluster 
positions with the APM galaxy survey, and for each cluster 
selected all galaxies falling within a distance of 1.2 h^^ Mpc 
from the cluster center. Since this is a larger radius than 
used in the cluster identification, some clusters near to the 
survey boundaries do not have complete data over the full 
circle. We found that 54 of the APM clusters are affected, 
and have simply rejected them from the sample, leaving 903 
clusters which we use in our analysis. 



3 PROJECTED CLUSTER SHAPES 
3.1 Basic Methods 

In order to estimate the APM cluster shapes we use the mo- 
ments of inertia method (cf. Carter & Metcalfe 1980; Plionis, 
Barrow & Frenk 1991) The gala^xy equatorial positions are 
transformed into an equal area coordinate system, centered 

on the cluster center, using: x = {Rag — Had) x coa{5ci) 
and y = 5g — Sd, where subscripts g and cl refer to galax- 
ies and the cluster, respectively. We then evaluate the mo- 
ments: hi = J2 '^ii'^i - ^i)' -^22 = J2 ^ii'^i ~ 2/<)' 

/i2 = hi = — WiXiUi, with Wi the statistical weight 
of each point. Note that because the inertia tensor is sym- 
metric wc have /12 = /21. Diagonalizing the inertia tensor 



det(7ij 



A'Ma) 



■ (M2 is 2 X 2 unit matrix.) 



(1) 



we obtain the eigenvalues Ai, A2, from which we define the 
ellipticity of the configuration under study by: e = 1 — A2/A1, 
with Al > A2. The corresponding eigenvectors provide the 
direction of the principal axis. 

This basic shape estimation method is applied using two 
alternative methods: 

(1) Discrete case [w = 1): In which we use the individual 
galaxies to determine the cluster shape. At first all galaxies 
within a small radius from the cluster center (~ 0.1 
Mpc) are used to define the initial value of the cluster shape 
parameters, then the next nearest galaxy is added consecu- 
tively to the initial group and the shape is recalculated until 
we include all galaxies within a limiting radius of our choice 
(usually ~ 0.75 hr^ Mpc). 

(2) Smooth case (w — 5): In which we use the smoothed 
galaxy distribution on a grid Ngr x Ngr, where the surface 
density of the j*^ cell is: 

J2iPj{Xi)F{Xi - Xgr) 



pjiXgr) — 



/ F{Xi - Xgr)A'^X 



(2) 



The sum is over the distribution of galaxies at positions Xi, 
with the Gaussian kernel being: 

{xi X gr^ 



F(x, 



2ttRs 



■ oxp 



27?? 



(3) 



All cells having a density above a chosen threshold are used 
to define the cluster shape, with cell-weights corresponding 
to the density fiuctuation: 

„„ _ Pjia^gr) - (p) 



iP) 



(4) 



where (p) is the mean projected APM galaxy density. Note 
that the apparent magnitude limit means that the APM 
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Figure 1. Distribution of maximum APM cluster radius used to 
define the cluster shape parameters for Rg = 0.24 (thick line) 
and Rs = 0.36 (hatched). Only the 405 APM clusters with no 
significant substructure are used (see section 3.4) 



clusters at large distances contain only the few most lumi- 
nous cluster galaxies. Therefore the smoothing radius Rsm 
must increase with cluster distance to obtain a continu- 
ous density field free of discreteness effects. We used sets 
of Monte Carlo cluster simulations with the same richness 
and depth distributions as the APM clusters to chose Ram, 
as a function of distance, so that it minimises discreteness 
effects and optimises the performance of our shape measur- 
ing algorithm. We find a roughly linear relation, fitted by: 
Rsm — 2.1 X 10~*r + 0.023, where r is the cluster distance. 

To estimate the cluster shape parameters we use cells 
above a threshold, defined as the mean (5) value of all grid- 
cells within a chosen radius (measured in Mpc), Rg, 
from the cluster center. After testing a range of values on 
mock clusters with the characteristics of the APM catalogue, 
we concluded that a value of Rg — 0.36 Mpc is a good 
compromise: too small a radius leads to inadequate sam- 
pling of the cluster region within ~ 1 Mpc; too large a 
radius extends the analysis to the very end of the sampled 
cluster region and thus artificially sphericalizes the clusters 
while it also increases the contamination effects by including 
relatively more background galaxies. Note that the highest 
density peak does not necessarily coincide with the listed 
APM cluster center. We have therefore measured moments 
about the position of the highest density peak, found within 
0.5 Mpc of the nominal centre. In only a few cases is 
there any significant difference between these two cluster 
centers. In any case, using either does not change the re- 
sults of our analysis, although on an individual basis it can 
have a marked difference to the derived cluster ellipticity, 
especially if the two center distance is large. 

In figure 1 we present the frequency distribution of the 
maximum cluster radius; ie. the distance between the cluster 
center and the most distant grid-cell that is used in the shape 
determination procedure, for two values of Rg. It is evident 
that there is a range of radii samples, depending slightly also 
on the value of Rg, the most common of which ~ 0.65 
Mpc. 



3.2 Testing the Performance of our Method 

The two procedures to determine cluster shape parameters 
may be biased by the unavoidable presence of foreground 
and background galaxies projected on the clusters, and may 
also have systematic biases due to the methods themselves. 
In this section we investigate the robustness of the two meth- 
ods in recovering the true projected cluster shape in the 
presence of discreteness effects and the galaxy background. 

To this end we generate a large number of mock clusters, 
resembling in appearance dynamically relaxed structures; 
ie., having no substructure and King- like surface brightness 
profiles: 



E(r) oc 



l+(^ 



(5) 



where Vc is the cluster core radius and a is the slope param- 
eter. Different values of the slope and core radius have been 
found in different studies, spanning the range 0.6^ 1 
and 0.1< rc< 0.25 h'^ Mpc (cf Bachall & Lubin 1993; Gi- 
rardi et al. 1995; Girardi et al. 1998). The latter study, using 
the ENACS optical sample of Abell clusters, found a median 



value of a ~ 0.65 



-t-0.05 



We create our mock clusters by ran- 



domly generating galaxies having a — 0.65, = 0.18 and 
input ellipticities of our choice. We have tested that small 
variations of these parameters do not alter significantly our 
results. 

The expected global background at each cluster can be 
estimated by: 



Nbac 



{p[z))z dz 



(6) 



where Zmax is the maximum redshift of galaxies in the APM 
catalogue (~ 0.3), d^li is the solid angle covered by the clus- 
ter, given by dQi = 27r(l — cosOi) for a cluster with angular 
radius 9i, and {p{z)) is the mean APM galaxy density at 
redshift z, obtained by integrating the APM galaxy lumi- 
nosity function, ${M,z), allowing for evolution (Maddox, 
Efstathiou & Sutherland 1996): 



${M,z)dM 



(7) 



with Mmin{z) the minimum absolute magnitude that a 
galaxy at a redshift z can have and still be included in the 
APM catalogue, limited in apparent magnitude by mum ~ 
20.5, ie; Mmin{z) = mum — 42.38 — 5 log z -I- 5 log /i — 3z . 
The points in figure 2 show the number of APM galaxies 
counted within a radius of 1.2 Mpc from each cluster 
center, and the line shows the predicted global background, 
Nback- Both are decreasing functions of distance due to the 
apparent magnitude limit of the APM galaxy catalogue. 

For our Monte-Carlo simulations we generate a ran- 
dom galaxy background, with number density given from 
eq.(|^, in a circular area of radius equal to the semi-major 
axis of the ellipse. As an example we plot in figure 3 the 
smooth galaxy density distribution of two mock clusters 
with eth = 0.5 at distances r = 200 Mpc and 380 
Mpc, respectively and having the typical APM cluster rich- 
ness at that distance. As expected the random background 
tends to sphericalize the clusters. 

The question that we want to answer now is: "Given 
an input cluster ellipticity what is the most probable mea- 
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Figure 2. The APM cluster richness (number of galaxies within 
a 1.2 Mpc radius) and the expected number of background 
galaxies as a function of distance. 
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Figure 3. The density contours for two mock clusters using a 
King-like profile and input ellipticity e^^ = 0.5. 



sured ellipticity recovered by our methods?" We present a 
case study of the mock cluster at a distance of 200 Mpc, 
having a range of input projected ellipticities of width 0.1. 
We generate 100 Monte-Carlo realizations of the cluster for 
each input ellipticity. In the left panel of figure 4 we present 
the measured mean ellipticities, and their scatter for both 
methods in the absence of a background galaxy distribution. 
Both methods give similar results, with the first method re- 
covering (for ethii 0.1) exactly the input (correct) elliptici- 
ties. 

The right panel of figure 4 presents the results for the 
same mock cluster but now we include the random back- 
ground galaxy distribution, appropriate for the distance of 
the cluster. The first method breaks down and severely un- 
derestimates the input ellipticities for e > 0.1. It is evident 
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Figure 4. Performance of our shape determination method for 
a simulated cluster at 200 Mpc containing 429 galaxies, by 
comparing the input and measured ellipticities, using either w = 1 
or w = 6 weights. The left panel presents the case for no back- 
ground galaxy contamination. The right panel presents the more 
realistic situation in which we superimpose a random background 
galaxy distribution according to eq. (^6h. 



that only the w — S method performs relatively well in the 
presence of the galaxy background and from now on we will 
be using only this method. 

Performing many tests with mock clusters at different 
distances, we conclude that the grid method has a variable 
performance as a function of distance, tending to increas- 
ingly underestimate the ellipticity of elongated clusters as a 
function of distance (by about 0.1-0.15). For clusters with 
eth~ 0.2 — 0.25 the relation between input and recovered el- 
lipticity is non-monotonic, which means it is not possible to 
correct the measured ellipticities of APM clusters for the sys- 
tematic biases which are evident in the right panel of figure 
4. These numerical tests have served mostly to choose be- 
tween the two shape determination procedures, rather than 
to derive a robust ellipticity correction procedure. 



3.3 Correcting systematic ellipticity biases 

Since cluster ellipticities cannot be corrected on an indi- 
vidual basis we must apply a statistical correction to the 
ellipticity distribution to deal with all the above mentioned 
systematic effects. In order to do this, we need to answer 
a slightly different question from the one we posed previ- 
ously. The relevant question is: "What is the distribution of 
input (correct) ellipticities from which a measured elliptic- 
ity can be obtained?" For each APM cluster with measured 
ellipticity, eoba, we generate 50 mock Monte-Carlo clusters 
for each ellipticity bin of width 0.05, spanning the whole 
range eth G (0, 1). These mock clusters, 1000 in total, have 
the same number of cluster and background galaxies and 
are placed at the same distance as the one observed. We 
measure the ellipticity of the mock clusters, emocfc, and de- 
rive its distribution function per each input (eth) ellipticity 
bin. We then measure, assuming Gaussian statistics, how 
many standard deviations the original APM cluster mea- 
sured ellipticity (tobs) differs from the derived mean mock 
value, (emock), and assign to the corresponding input ellip- 
ticity (eth) bin the resulting probability. Therefore, for each 
input (eth) ellipticity bin i we estimate the probability of 
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Figure 5. The probability functions for some characteristic mea- 
sured eUipticities (tobs = 0.2,0.5 and 0.7) of three APM clusters 
containing ~ 140 galaxies and at a distance of ^ 320 Mpc. 



being the correct ellipticity from which the APM measured 
one could have resulted; ie, 

P'(e„0 = 1 - P{z,) (8) 

where Zi = \eohB — {^\nock)\l '^i- If for example, for some bin 
i, the value of |eobs — (cmocfc)! = 4(Ti then Zi = 4, P{zi) ~ 1 
and thus the probability of ej^ being the correct APM clus- 
ter ellipticity is ~ 0. Doing so for all bins we derive the 
full probability distribution function of input (correct) eUip- 
ticities, P{eth), from which our measured APM cluster el- 
lipticity could have resulted. For each APM cluster we have 
therefore generated in total 1000 mock clusters (50 per eth 
bin) which provide an estimate P{eth) for that cluster. The 
final step in estimating the corrected cluster ellipticity dis- 
tribution is to add the contribution from each cluster by ran- 
domly sampling its P(eth) A*' times (we have used N = 20). 
We stress that the correction procedure is applied separately 
to each individual APM cluster which means that the cor- 
rections take into account the variations in performance of 
our shape determination method as a function of different 
cluster richness and distance. 

As an illustrative example we plot in figure 5 the P{^th) 
distribution for three distant clusters (r ~ 320 Mpc) 
with eobs — 0.2, 0.5 and 0.7 respectively. For the eobs = 0.2 
case the P{(.th) is quite flat in the range 0^ ethi$ 0.2, hav- 
ing a long tail up to eth — 0.6. The eobs ~ 0.5 case has a 
P{eth) distribution that peaks at ~ 0.65 but also has a sig- 
nificant contribution from lower eUipticities. Finally for the 
^obs = 0.7 case the most probable value of eth is ~ 0.9. These 
facts serve to show that in order to produce the correct ap- 
parent cluster ellipticity distribution it is essential to use 
a statistical procedure and sample each P{eth) distribution 
adequately. 



3.4 Substructure of APM clusters 

The correction procedure of the ellipticity distribution is 
based on the assumption that the clusters under study have 
a smooth King-like density profile. Therefore we need to 
exclude from our sample clusters that exhibit evidence of 
significant substructure. The number of clusters with sub- 
structure, expected to be undergoing dynamical evolution, is 
an unsettled issue but of great importance since it provides 
information of the mean density of the universe (cf. Rich- 
stone et al. 1992; Dutta 1995; Buote 1998; Thomas et al. 
1998). In an = 1 universe, clusters continue to form even 
today and therefore one expects more substructure than in 
a low Q universe. Identifying real sub-clumps in clusters is a 
difficult problem in general, since one has to work either in 
two-dimensions, in which projection effects can significantly 
affect the visual structure of clusters, or in redshift space, 
where distortions due to knowledge of only the radial veloc- 
ity component can again distort the true pattern. Several 
studies indicate that at least ~ 30% - 50% (cf. West 1994 
and references therein; Jones & Forman 1999) of rich clusters 
have strong substructure in their gas or galaxy distribution 
within ~ 1 Mpc of the cluster center. 

Here, we present the main steps of our substructure 
identification method, the full details and results will be pre- 
sented in a forthcoming paper. We work on the smoothed 
density field, as described in section 3.1. For each overden- 
sity threshold, estimated for Rs = 0.12,0.24 and 0.36 
Mpc, we select all grid-cells with overdensities above the spe- 
cific threshold. We then connect all cells that have common 
borders to create multiple clumps. In all cases we accept 
only clumps that are within ~ 0.75 Mpc of the highest 
cluster peak; we have found this scale to be optimal for re- 
liable substructure identification by validating our methods 
and results using ROSAT data for a subsample of 22 clus- 
ters (Kolokotronis et al. 2000), which also corresponds to the 
counting radius used in the APM cluster finding algorithm 
(cf. Dahon et al. 1997). 

Investigation of the number and size of these clumps as 
a function of overdensity threshold, provides the following 
categorisation: 

• No substructure (69 clusters): Clusters with one clump 
in all overdensity levels. 

• Weak substructure (336 clusters): Multiple clumps only 
at the lowest overdensity level or at the highest two over- 
density levels but where the second in size clump is < 20% 
of the total cluster size (cf. Richstone et al. 1992). 

• Strong substructure (498 clusters): Multiple clumps 
where the second in size clump is ^ 20% of the total cluster 
size. 

We have investigated the robustness of our substructure 
characterisation procedure, as a function of different Rs, and 
found only small variations, consisting mainly of a move- 
ment of APM clusters between the "no" and "weak" sub- 
structure categories. We have also verified that due to the 
random galaxy background, coupled with discreteness ef- 
fects, it is common to find "weak" substructure even in our 
mock clusters which by construction have no substructure. 
Therefore we exclude from our shape determination analysis 
only those APM clusters that were found to have "strong" 
substructure. 
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4 TRUE CLUSTER SHAPES 

In order to find the intrinsic ellipticity distribution assuming 
tiiat clusters arc all oblate or prolate spheroids, we use a 
standard method based on the kernel estimator. 



4.1 Kernel estimator 

General reviews of kernel estimators are given by Silver- 
man (1986) and Scott (1992) but the applications to the 
astronomical data are given by Vio et al. 1994, Tremblay & 
Mcrritt (1995) and Ryden (1996). Here we review the ba- 
sic steps of the Kernel method, following the notation of 
Ryden (1996). For each estimated ellipticity e, we estimate 
the axial ratio, q = (-j^)^^^. Given the sample of axis ra- 
tios qi,q2, ■■■■,qN for A'^ clusters, the kernel estimate of the 
frequency distribution is defined as: 



(9) 
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where K{t) is the kernel function, defined so that 

K{t)dt = 1 , (10) 



/ 



and h is the "kernel width" which determines the balance 
between smoothing and noise in the estimated distribution. 
In general the value of h is chosen so that the expected value 
of the integrated mean square error between the true, f{q), 
and estimated, f{q), distributions, ^fK{x) — /(a;)] dx, 
is minimised (cf. Vio et al. 1994; Tremblay & Merritt 1995). 
In this work we estimate the h using a very common ap- 
proach presented by Silverman (1986), Vio et al. (1994) and 
Ryden (1996) in which: 



h = Q.9AN~'^f^ 



(11) 



where A'^ is the number of the clusters and A = 
min(a", (54/1.34), with Q4, the interquartile range. There are 
three common choices for the kernel function K{t) which 
have quadratic, quartic and Gaussian forms (cf. Tremblay 
& Merritt 1995). Many studies have shown that the choice 
of a kernel function does not in general affect the estimates, 
and they differ trivially in their asymptotic efficiencies. We 
have chosen a Gaussian kernel: 
1 



2tt 



(12) 



In order to obtain physically acceptable results with 
f(q) =0 for q < Q and q > 1, we apply reflective boundary 
conditions which means that the Gaussian kernel is replaced 
with: 



2 - q - Qi 



h 



(13) 



Kiq,q.M-K{l^).K{i±-),K{ 

This also ensures the correct normalization, J^^ f{q)dq = 1. 
For a discussion of reflective boundary conditions see Silver- 
man (1986) and Ryden (1996). The crosses in figure 6 show 
the projected axial ratio distributions with the Poisson Icr 
error bars and the solid lines show the kernel estimate / with 
width h = 0.075. In the top panel we present our results for 
the uncorrected distribution of all 903 APM clusters; which 
can crudely be fitted by a Gaussian with (g) ~ 0.65 and 



Figure 6. In tiic top panel the crosses show the uncorrected 
distribution of apparent axial ratios for the complete sample of 
903 APM clusters with the solid line being the smooth distribu- 
tion from the nonparametric kernel estimator. The open circles 
correspond to the distribution of the 405 APM clusters with no 
significant substructure. In the lower panel we present the cor- 
rected axial ratio distribution for the later clusters. 



standard deviation ~ 0.15. In the bottom panel we present 
the corrected distribution using the 405 APM clusters free 
of significant subclustering. It is obvious that the two dis- 
tributions are significantly different, with the peak of the 
corrected distribution having moved to lower q's but with 
an extended contribution of apparently quasi-spherical ob- 
jects. 



4.2 Inversion method 

The relation between the apparent and intrinsic axial ra- 
tios, is described by a set of integral equations first inves- 
tigated by Hubble (1926). These are based on the assump- 
tions that the orientations arc random with respect to the 
line of sight, and that the intrinsic shapes can be approx- 
imated by either oblate or prolate spheroids. There is no 
physical justification for the restriction to oblate or prolate 
but it greatly simplifies the inversion problem. Furthermore, 
if the intrinsic shape of clusters is triaxial or a mixture of the 
two spheroidal populations then there is no unique inversion 
(PBF). Writing the intrinsic axial ratios as /3 and the esti- 
mated distribution function as No{fi) for oblate spheroids, 
and Np (/3) for prolate spheroids then the corresponding dis- 
tribution of apparent axial ratios is given for the oblate case 
by: 



/(?) 



Jo 



N.{l3)dl3 



(1 _g2)l/2(^2 _^2)l/2 

and for the prolate case by: 

^^"^ (l-g2)l/2(g2_ ^2)1/2 



(14) 



(15) 
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Figure 7. The distribution of intrinsic APM cluster axial ra- 
tios and its la range, for the corrected and uncorrected samples, 
assuming that clusters are either oblate or prolate spheroids. 



Inverting equations (eq.|l^ and (eq.^^ gives us the distri- 
bution of real axial ratios as a function of the measured 
distribution: 



and 



, dg (/32-g2)l/2 



7r/3 



(16) 



(17) 



with /(O) = 0. In order for iVp(/3) and No(P) to be phys- 
ically meaningful they should be positive for all /3's. Fol- 
lowing Ryden (1996), we numerically integrate eq.(|l^) and 
eq.(|l7|) allowing Np{P) and No{f3) to take any value. If the 
inverted distribution of axial ratios has significantly negative 
values, a fact which is unphysical, then this can be viewed 
as a strong indication that the particular spheroidal model 
is unacceptable. 

In figure 7 we present the uncorrected and corrected 
intrinsic axial ratio distributions. The uncorrected distribu- 
tion for both spheroidal models takes negative values; over 
the range fJ^ 0.7 for the oblate case, and 0.8 for the 
prolate case. Using the corrected apparent axial-ratio dis- 
tributions f{q), the oblate model produces negative values 
of N for /?> 0.5 and /?< 0.2, but the prolate one provides a 
distribution of intrinsic axial ratios that is positive over the 
whole P range. 

This suggests that the APM cluster shapes are better 
represented by that of prolate spheroids rather than oblate, 
which is in agreement with PBF and Cooray (1999). How- 
ever, it is probably not realistic to assume a population of 
pure oblate or prolate spheroids but rather of triaxial ellip- 
soids, in which case the inversion procedure is not unique 



(see PBF). Nevertheless, our results strongly suggest that 
cluster prolateness should be a dominant feature. 



5 CONCLUSIONS 

We have measured the projected ellipticities of all APM clus- 
ters using moments either of the individual galaxy distribu- 
tion or of the smoothed galaxy distribution above some over- 
density threshold. We have performed large sets of Monte 
Carlo simulations in order to test the statistical robustness 
of the two procedures, and conclude that the first method 
is strongly affected by the presence of background galaxies 
whereas the second method better recovers the underline 
true cluster projected ellipticity. 

We devised a statistical Monte-Carlo procedure to cor- 
rect the distribution of cluster ellipticities for the system- 
atic errors introduced by discreteness effects, background 
galaxies and the method itself. The procedure involves es- 
timating first the probability distribution of the true pro- 
jected ellipticity for each APM cluster, and then random 
sampling this probability distribution to estimate the clus- 
ter's contribution to the corrected overall projected cluster 
ellipticity distribution. This method works well only for clus- 
ters that appear relaxed, with no significant substructure. 
'Therefore we have excluded, from our final cluster sample 
which contains 405 clusters, all the APM clusters that show 
evidence of significant substructure. Prior to the exclusion 
of these clusters the uncorrected axial ratio distribution of 
the whole APM cluster sample can be crudely approximated 
by a Gaussian with a mean of ~ 0.65 and a standard de- 
viation of ~ 0.15. The corrected apparent axial- ratio distri- 
bution is significantly different showing a bump at g ~ 0.46 
and having a significant contribution from apparently quasi- 
spherical systems. 

Using the nonparametric kernel procedure we obtain a 
smooth estimate of the apparent APM cluster axial-ratio 
distribution. We assume that the APM clusters are a ho- 
mogeneous population of either oblate or prolate spheroids 
and numerically invert the apparent distribution to obtain 
the intrinsic distribution. The most acceptable model is pro- 
vided by that of prolate spheroids. This result supports the 
view by which clusters form by accretion of smaller units 
along the large-scale structure (filament) in which they are 
embedded (cf. West 1994; West, Jones & Forman 1995). 
Such an accretion process would happen preferentially along 
the cluster major axis, which is typically aligned with the 
nearest cluster neighbour (cf. Bingelli 1982; Plionis 1994 and 
references therein). 

Since cluster shapes and substructure are sensitive cos- 
mological probes (Evrard et al. 1993), we plan to investi- 
gate these issues further using APM clusters and compare 
our results with theoretical expectations, provided by high- 
resolution N-body simulations of different cosmological mod- 
els (cf. Thomas et al.l998). 
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